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ABSTRACT 

The influence of inelastic material models on computed stress-strain 
states, and therefore predicted lives, was studied for thermomechanical ly 
loaded structures. Nonlinear structural analyses were performed on a 
fatigue specimen which had been subjected to thermal cycling in fluidized 
beds and on a mechanically load-cycled benchmark notch specimen. Four in- 
cremental plasticity-creep models (isotropic, kinematic, combined isotropic- 
-kinematic, combined plus transient creep) were exercised using the MARC 
program. Of the plasticity models, kinematic hardening gave results most 
consistent with experimental observations. Life predictions using the com- 
puted strain histories at the critical location with a Strainrange Parti- 
tioning approach considerably overpredicted the crack initiation life of the 
thermal fatigue specimen. 

INTRODUCTION 

Hot section components of aircraft gas turbine engines, such as combus- 
tor liners and turbine blades and vanes, are subject to progressive creep- 
fatigue damage resulting from cyclic thermomechanical loading under extreme 
gas pressure and temperature environments. A Strainrange Partitioning ap- 
proach (ref. 1) to assess the durability of these components has been under 
development at the NASA Lewis Research Center. In order to apply this or 
similar methods, it is first necessary to determine the stress-strain- 
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temperature history of the part at the critical location where cracks will 
initiate. 

As part of the life prediction studies at Lewis, wedge specimens have 
been thermally cycled in fluidized beds as described in reference 2. In 
these tests, two fluidized beds were used to rapidly heat and cool prismatic 
bar specimens of single or double edge wedge cross-section. Nonlinear 
structural analyses were performed for these specimens using the MARC pro- 
gram (ref. 3); the results are reported in references 4 and 5. These non- 
linear analyses were for specimens of several alloys and used a combined 
isotropic-kinematic hardening model in MARC in conjunction with monotonic 
stress-strain properties taken from the literature. 

Finite-element nonlinear analysis methods are becoming of increasing 
interest for computing the cyclic stress-strain response of hot section com- 
ponents (refs. 6 to 10). A major disadvantage of these methods, excessive 
computing costs, is being alleviated by advances in computer technology. 
Another deficiency is that current nonlinear analysis computer codes utilize 
classical constitutive material models whose accuracies vary with the type 
of material and the cyclic conditions involved. Furthermore, these class- 
ical models simplify the analyses by uncoupling time independent (plas- 
ticity) and time dependent (creep) effects, neglecting strain rate effects 
on plastic flow, and defining specific yield surfaces. The NASA Lewis 
Research Center has instituted programs to develop constitutive models which 
would more realistically represent the inelastic material behavior and be 
computationally practical for finite-element structural analysis. To verify 
the nonlinear structural analysis methodologies, Lewis is also sponsoring 
controlled cyclic experiments to provide strain data for benchmark notch 
specimens (ref. 11). 
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In this study existing constitutive models in the MARC computer program 
were exercised in inelastic analyses of an IN 100 wedge specimen subjected 
to thermal cycling and an Inconel 718 benchmark notch specimen subjected to 
mechanical load cycling. The objective of the study was to evaluate the 
effects on calculated hysteretic response, and therefore predicted life, of 
different inelastic constitutive models available in nonlinear analysis com- 
puter codes. 

Three dimensional elastic and nonlinear structural analyses were per- 
formed on a thermally cycled double-edge wedge specimen. The nonlinear an- 
alyses were conducted using isotropic, kinematic and combined isotropic- 
kinematic hardening models and a combined hardening model in conjunction 
with a strain hardening creep law to account for cyclic time-dependent ef- 
fects. Strain histories computed at the critical location from the dif- 
ferent constitutive models were used in conjunction with the Strainrange 
Partitioning method to compare predicted lives against the observed crack 
initiation life. Two dimensional nonlinear analyses were performed for a 
mechanically load-cycled benchmark notch specimen; computed strain histories 
at the notch root using various material models were compared against 
measured notch strains. 

PROBLEM DESCRIPTION 

The primary structure considered in this study was an IN 100 alloy 
double-edge wedge specimen as illustrated in figure 1. Cracking was ob- 
served at the 1/4 span position on the leading edge after 38 cycles of test- 
ing in the fluidized bed facility (ref. 2). 

The physical properties of the cast IN 100 alloy are presented in 
table I. Mean thermal coefficient of expansion data were converted to in- 
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stantaneous coefficients of thermal expansion for input into the MARC pro- 
gram. The modulus of elasticity was determined from monotonic stress-strain 
tests of tensile specimens. Cyclic stress-strain curves were obtained using 
the single specimen incremental step procedure and equipment described in 
reference 12. A typical cyclic stress-strain curve, with the loci of the 
curve tips represented by an exponential equation, is illustrated in figure 
2. Also shown for comparison in figure 2 is a monotonic stress-strain curve 
represented by an exponential equation. Short-time cyclic creep tests were 
conducted on IN 100 specimens using the procedures and facilities described 
in reference 13. Preprocessor programs expressed both the cyclic stress- 
strain and creep data as functional relations in exponential form. These 
equations were incorporated into MARC by means of user subroutines. The 
constants of the cyclic and monotonic stress-strain equations are given in 
table II for various temperatures. In table III the constants of the cyclic 
creep equations are given for various temperatures. 

The specimen was thermally cycled in fluidized beds maintained at 316° 
and 1088° C with an immersion time of 3 minutes in each bed. Transient tem- 
perature loading on the specimens was determined from thermocouple data as 
described in reference 2. Curve fits of thermocouple data along the mid- 
chord at the midspan at various increments after immersion into the flui- 
dized beds are presented in figure 3. The temperature gradient through the 
thickness of the wedge was assumed to be negligible. Another set of thermo- 
couple data was taken with thermocouples mounted along the leading edge over 
half the span to obtain the longitudinal (along the span of the specimen) 
temperature gradient for the different time increments. 

Supplemental analyses to evaluate the constitutive material models were 
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also performed for a benchmark notch specimen of Inconel 718 alloy which was 
load cycled at a frequency of 0.167 Hz and a temperature maintained at 
649° C. The material properties given in reference 11 were correlated in 
the same way as the IN 100 alloy properties. 

ANALYTICAL PROCEDURE 

Stress and total-plastic-creep strain distributions in the wedge speci- 
mens were calculated from the MARC nonlinear, finite-element computer pro- 
gram. Computations were performed for 34 time increments (17 heating, 

17 cooling) into which the thermal cycle was subdivided, as shown in fig- 
ure 3. The analyses were terminated when stable stress-strain hysteresis 
loops were obtained or after three cycles if the hysteresis loops remained 
unstable. 

Plasticity computations were based on incremental plasticity theory us- 
ing the von Mises yield criterion and normality flow rule. The yield sur- 
face under reversed loading was determined from the stress-strain properties 
and the selected hardening model. Three hardening models available in MARC 
(isotropic, kinematic and combined isotropic-kinematic) were selected for 
evaluation. Monotonic stress-strain properties were used in conjunction 
with the isotropic and combined models because of their initial insta- 
bility. Saturated cyclic stress-strain properties were used for the stable 
kinematic model. A bilinear representation of the cyclic stress-strain 
curve, as shown in figure 2, was applied to the kinematic hardening model. 
The slope of the kinematic model was determined from energy considerations 
so that the strain energy, as indicated by the enclosed area, would be iden- 
tical with that of the actual cyclic stress-strain curve. Creep effects 
during the cycle were considered for one case involving the combined model 
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by imposing four 30 second hold times during heating and two 6 second hold 
times at the start of the cooling part of the cycle. These intervals were 
selected because the combination of temperatures and stresses indicated a 
possibility of the occurrence of significant creep at these times in the 
thermal transient. The creep computations utilized the cyclic creep data in 
conjunction with a strain-hardening rule. A subroutine which was inserted 
into the MA.RC program in the form of yield strengths and work hardening 
slopes as functions of temperature, was used to determine the stress-strain 
properties for the local temperatures at the Gaussian integration points. 
Similarly the creep properties and laws were coded into another user sub- 
routine which was used to obtain the creep strains at the integration points. 

A preprocessor program converted the thermal loading data from the wedge 
specimen into the form of sixth-order polynomial equations. A subroutine, 
which was inserted into MARC, interpolated from these equations for the 
local temperatures at the Gaussian integration points. 

The finite element model for the wedge specimen is illustrated in fig- 
ure 4. Because of symmetry only one-fourth of the specimen needed to be 
modelled; this model was bounded by the surface and intersecting midchord 
and midspan planes of symmetry. The element used was a 20 node, isopara- 
metric, three dimensional block with 8 corner nodes and 12 edge midpoint 
nodes. This element had 27 Gaussian integration points. The model con- 
sisted of 36 elements with a total of 315 nodes and 778 unsuppressed degrees 
of freedom. 

All nodes initially on the midspan and midchord faces of the model were 
constrained to lie on the midspan and midchord planes respectively. In ad- 
dition, one node at the leading edge was constrained chordwise (leading to 
trailing edge) in order to prevent rigid body motion in that direction. 
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The analytical procedure used for the benchmark notch specimen was ba- 
sically the same as for the wedge specimen. Each cycle was subdivided into 
30 load steps. One fourth of the specimen was modelled as shown in figure 5 
using 592 plane strain, triangular elements with a total of 335 nodes. 

RESULTS AND DISCUSSION 

The critical location for crack initiation in the thermally cycled 
double-edge wedge specimen is at the leading edge at a quarter of the speci- 
men span from either end. Results of both elastic and inelastic structural 
analyses determined that the critical location based on the region of the 
finite element model with the largest total strain range during the cycle 
was coincident with the observed crack initiation site. In the following 
discussion, the stress-strain results for the critical location were actu- 
ally computed at the closest Gaussian integration point which was 0.056 cen- 
timeter from the surface at the quarter span. 

The stress-total strain solutions at the critical location from the MARC 
elastic and nonlinear analyses of the wedge specimen are shown in figure 6. 
All stresses and strains in this figure were effective or equivalent values 
which were originally computed as positive numbers. However, in order to 
construct stress-strain hysteresis loops for life prediction purposes, the 
stresses and total strains were assigned positive or negative signs depend- 
ing on the signs of the highest magnitude principal stresses or strains. 
Nonlinear stress-strain hysteresis loops are presented for the second cycle 
of the analyses. During heating the metal temperature at the critical loca- 
tion increased from 343° C at the start of the cycle to 1077° C at the end 
of heating. In all analytical cases, the minimum total strain occurred 
after 30 seconds of heating when the temperature at the critical location 
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was 888° C and the maximum total strain after 9 seconds of cooling or a 
total elapsed time of 189 seconds when the temperature was 749° C. 

Predicted stress-strain hysteresis loops from the elastic analysis and 
the nonlinear analyses using combined and kinematic hardening models are 
compared in figure 6(a). These results indicate that the total strain range 
was not appreciably affected by the choice of constitutive model or type of 
stress-strain data and that an elastic analysis was adequate for the compu- 
tation of the total strain range. The major differences between the elastic 
and nonlinear hysteresis loops were in the stress levels, which shifted in 
the tensile direction under inelastic straining with the largest peak and 
mean stresses obtained with the combined hardening model. A measure of the 
strain energy or plastic work is the area of the hysteresis loop. The 
widest hysteresis loop and, therefore, the most plastic work is shown by the 
kinematic hardening model in figure 6. There was no further plastic strain- 
ing or work during or after the second cycle using the combined hardening 
model and, therefore the area and shape of the combined and elastic hyste- 
resis loops in figure 6 remained about the same. 

The nonlinear analysis using the isotropic hardening model gave essenti- 
ally the same stress-strain solutions as were obtained with the combined 
model in figure 6(a) due to the use of the same monotonic stress-strain 
properties and the absence of plastic strain reversal during cycling. 
Therefore, the discussion of results for the combined hardening model is 
also applicable to isotropic hardening and the latter will not be discussed 
separately. 

Figure 6(b) compares the stress-strain hysteresis loops from the non- 
linear analyses using the combined hardening models with and without creep. 
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Inclusion of creep effects during the thermal transients had only a small 
effect on the peak and mean stresses with combined hardening, but resulted 
in substantially more strain energy per cycle as represented by the enclosed 
areas of the stress-strain hysteresis loops. Although the hysteresis loops 
for the combined-creep and kinematic models in figure 6 are shown as closed, 
there was some inelastic strain ratchetting which was relatively minor and 
therefore ignored in plotting the loops. 

Stabilization of the stress-strain solution using the combined hardening 
model is shown in figure 7(a) where it is seen that there was no further 
plastic flow after the first 60 seconds of heating; this is an impossibility 
since the specimen could not fail in 38 cycles without undergoing sub- 
stantial plastic strain cycling. In contrast the kinematic hardening re- 
sults in figure 7(a) exhibit plastic strain reversal and ratchetting with a 
relatively constant plastic strain range per cycle. Figure 7(b) shows the 
inelastic strain response for the combined-creep case. Accounting for tran- 
sient creep effects resulted in creep strain ratchetting on every cycle and 
smaller plastic strain changes with the combined hardening model. Only 
slight changes in the maximum equivalent creep strain were obtained with 
further cycling. However, the minimum equivalent creep strain increased, 
and therefore the creep strain range decreased, although at diminishing 
rates during cycling. 

The computed strain histories at the critical location were used to pre- 
dict crack initiation life based on the Strainrange Partitioning Life Pre- 
diction Method. The material life relationships for this method are defined 
in reference 14 for cast IN 100 alloy from isothermal fatigue and creep rup- 
ture tests. For these analyses the response from the kinematic model con- 
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tained only pp(tensi1e plasticity reversed by compressive plasticity) and 
from the combined-creep model was conservatively assumed to contain cc(ten- 
sile creep reversed by compressive creep) damage cycles. Crack initiation 
lives of approximately 1400 cycles were predicted in both cases compared to 
the observed life of 38 cycles. The overpredictions in life are not neces- 
sarily proof of the inadequacy of the structural analysis method since there 
is evidence that thermal cycling produces damage at a faster rate than com- 
parable isothermal, strain-controlled test data used in the life prediction 
method. 

In figure 8 analytical results using both combined and kinematic harden- 
ing models are compared against the experimental load-notch strain cycle 
from the benchmark notch test. Creep was not a significant factor under the 
continuous cycling, isothermal conditions of this test. The experimental 
results demonstrated that a stable load strain response occurred on the 
first cycle with only minor strain changes due to subsequent cycling. A 
plasticity analysis using the combined hardening model did not accurately 
represent the experimental results; it predicted, after initial loading, an 
elastic response with further cycling (fig. 8(a)). Another plasticity anal- 
ysis using the kinematic hardening model demonstrated good agreement with 
the experimental results. Kinematic hardening predicted ratchetting between 
the first and second cycles and a stable notch strain cyclic response there- 
after (fig 8(b)); except for slightly overpredicting the ratchetting, these 
results are consistent with the experimental notch cyclic response. 

SUMMARY OF RESULTS 

The results of the evaluation of inelastic constitutive models available 
in nonlinear, structural analysis computer programs can be summarized as 
follows: 
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1. Of the plasticity hardening models which were evaluated, the kin- 
ematic model gave a predicted stress-strain response most consistent with 
experimental observations. The combined (as well as the isotropic) model 
predicted elastic response during cycling which obviously did not agree with 
experimental results from both the thermal fatigue wedge and benchmark notch 
specimen tests. Creep effects were shown to be significant during thermal 
transients and failure to take them into account can affect the predicted 
stress-strain response. 

2. Of the structural analysis parameters used in low-cycle fatigue dam- 
age models only the total strain range was relatively insensitive to the 
choice of inelastic constitutive model. Other parameters such as inelastic 
strain range, mean stress, and inelastic work were significantly affected by 
the constitutive model. The elastically computed maximum total strain range 
agreed well with that computed from the inelastic analyses. The elastic 
analysis was also able to determine the critical location for crack initi- 
ation and the cycle times when the total strain was maximum or minimum. 

3. The life prediction analyses based on the structural analysis results 
using the kinematic and combined-creep models in conjunction with iso- 
thermal, strain controlled fatigue test data overpredicted the observed 
crack initiation life of the thermally-cycled wedge specimen. 
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TABLE 1 


IH 100 ALLOY PHYSICAL PROPERTIES 



Modulus of Elasticity, 1 Mean Coefficient of Thermal Expansion, I 

»C j 

MH/n* 1 

m/m/®C 

316 i 

193X103 i 

13.1X10"* 

371 1 

190 1 

13.3 

627 1 

186 i 

13.5 

682 1 

183 1 

13.7 

538 1 

179 1 

13.9 

593 1 

176 1 

16.0 

669 1 

172 1 

16.6 

706 1 

168 1 

16.6 

760 1 

163 1 

16.9 

816 1 

157 1 

15.6 

871 1 

152 1 

15.8 

927 1 

165 1 

16.6 

982 1 

139 1 

16.7 

1038 1 

133 1 

17.5 

1093 1 

127 1 

18.2 


TABLE II. IMIOO ALLOY STRESS-STRAIH PROPERTIES 


Temperature, 

oc 

1 

1 i ‘ 

Cyclic, 

cf=Kle,y 

.1)" 

1 

1 

Monotonic, 
'X a=C(c^/.l)"> 

1 K 

j 

n 

j 

C 

m 

316 

1 1005 

i 

.066 

1 

731 

.078 

627 

1 966 

1 

.066 

1 

731 

.078 

538 

1 869 

1 

.086 

1 

731 

.078 

669 

1 777 

1 

.113 

1 

731 

.078 

760 

1 665 

1 

.167 

1 

731 

.078 

871 

1 528 

1 

.187 

1 

676 

.078 

982 

1 361 

1 

.236 

1 

255 

.166 

1093 

1 157 

1 

.297 

1 

173 

.166 


• Locus of cyclic curve tips (fig. 2) 

> Stress (<t) in MPa, plastic strain (e^) in percent 
X Not applicable for 6^ less than 0.02 percent 


TABLE III. IH 100 AHOY CREEP PROPERTIES 


Temperature# i Creep rate, X/min., 

OC I ‘ 6 =A(ff/5.895)"'(t)" 

I _ 


■■■IH 

■■■■Qlll 

n 

m 

1 n 

760 

1 

1 

1 

1 

1 .00062 

1 

1 

1 

1 

I 

0.717 

1 -0.881 

871 

1 .00012 

1 

1.709 

1 -0.736 

982 

1 .00010 

I 

2.172 

1 -0.656 

1093 

1 .00058 

I 

2.103 

1 -0.636 


• Stress (ff) in MPa, time (t) in minutes 


102 



















Trailing- 
edge radius, 


a 102 ■ 


1.635 


a 317 


4 '^ 

^ K3.17H '''■ 


a 635 


aawH 


y? '^Leading-edge 
radius, 0,066 


la 16 


Figure 1. - Double-edge wedge. (All linear dimensions 
in centimeters. I 
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Figure 2. - Representation of stress-strain curves. 
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Figure 3. - Temperature of midchord it midspan at various times after Immersion Into fluidized beds. 
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Figure 4. - Model and typical element used (or MARC . Benchmark notch specimen finite 

analysis with coordinate convention. element model (all dimensions in cml. 
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(a) Comparison of combined, kinematic models and elastic results. 


(b) Comparison of combined and combined-creep models. 


Figure 6. - Effect of Inelastic constitutive model on stress-strain response 
of wedge specimen (critical location). 
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(b) Combined-creep model. 


Figure 7. - Inelastic strain response at critical location o( wedge 
specimen. 




Figure 8. - Comparison of benchmark notch specimen experi- 
mental and analytical results. 
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